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ABSTRACT 

Context. While most sources above 10 5 L Q have already formed an Ultracompact Hn region (UCHu), this is not necessarily the case for sources 
of lower luminosity. Characterizing sources in the transition phase, i.e., very luminous objects without any detectable free-free emission, is 
important for a general understanding of massive star formation. 

Aims. Characterizing one of the most luminous High-Mass Protostellar Objects (HMPO) that has not yet formed any detectable UCHu region. 
Methods. The region was observed with the Submillimeter Array in three different array configurations at 875 /im in the submm continuum 
and spectral line emission at sub-arcsecond resolution. 

Results. The 875 fim submm continuum emission has been resolved into at least two condensations. The previously believed driving source of 
one of the outflows, the infrared source IRS1, is ~ 0.9" offset from the main submm peak. The data do not allow to differentiate whether this 
offset is real, either caused by different sources or a shift of the photo-center due to scattering, or whether it is only due to poor astrometry 
of the infrared data. Over the entire 4 GHz bandwidth we detect an intermediate dense spectral line forest with 27 lines from 8 different 
species, isotopologues or vibrationally-torsionally excited states. Temperature estimates based on the CH3OH line series result in values of 
T(Peakl) ~ 150 ± 50 K and T(Peak2) ~ 80 ± 30 K for the two submm peak positions, respectively. The SiO(8-7) red- and blue-shifted line 
maps indicate the presence of two molecular outflows. In contrast, the vibrationally-torsionally excited CH3OH line exhibits a velocity gradient 
approximately perpendicular to one of the outflows. With a size of approximately 5000 AU and no Keplerian rotation signature, this structure 
does not resemble a genuine accretion disk but rather a larger-scale rotating toroid that may harbor a more common accretion disk at its so far 
unresolved center. 
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1. Introduction 

In contrast to known Ultracompact Hn regions (UCHns) 
with luminosities sometim es even exceeding 10 6 L o (e.g., 
Wood & Churchwelllll989l) . a close inspection of large sam- 
ples of younger High-Mass Protostellar Objects (HMPOs) 



(e.gjMolinari et al.ll996tlSridharan et al.ll2002l:lFaundez et al 
l2004t IBeltran et al.l 120061) shows that only very few sources 
at that evolutionary stage exceed 1O 5 L . The average lumi- 
nosity difference between HMPOs and UCHns is likely at- 
tributed to evolutionary differences, i.e., the still accreting 
HMPO will gain more mass eventually being energetic enough 
to produce detectable free-free emission at c m wavelengths 



and hence evolve in t o an U CHu region (e.g., Sridha ran et al 
2002; Beuther et al. 2007a). Observations indicate that be- 
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low 10 5 L the free-free emission of the forming star is still 
potentially quenched, either by the pressure of the infalling 
gas, high recombination rates due to hi gh accretion and infall 
rates, or by gravitational trapping (e.g. JWalmslevlll995l ; iKeto 
2003). However, since almost all massive star-forming regions 
above 10 s L are detected at cm wavelengths, these processes 
appear less feasible in that regime. Most interferometric in- 
vestig ations focused either on ve r y lum i nous UCHu regions 
(e.g.. ICesaroni et alj 11994 iKetol l2002t IBeltran et alj [2004; 



Hoare et al 



2007) or on younger but less luminous HM POs 



(e.g. [Shepherd et al. 1998; Zhang et alj|1998t ICesaroni et al 



1999; Beuther et alj|2002ai l2005ch . This transition phase of a 



source that has already a very luminous and massive accret- 
ing protostar, however, that is yet not capable to ionize enough 
gas to be detected as a Hypercompact Hn region (HC11), has 
not been studied well so far. Therefore, here we present new 
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Submillimeter Array (SMAQ) data observed toward the 10 L 
HMPO IRAS 23151+5912. 

The rationale to choose IRAS 23151+5912 was foremost 
based on two criteria: the source had to be very luminous 
(> 10 5 L o ) and it had to be very young, in an evolution- 
ary stage prior to forming any detectable ultracompact Hn re- 
gion. The source IRAS 23151+5912 comprised both features, 
hence it is one of the most luminous HMPOs in a pre-UCHn 
region phase. The source at a kinematic distance of 5.7 kpc 
is part of a large sample of 69 c andidate HMPOs init ially 
studied by Sri dharan et al. (120021) : Beuth er et alj I (l2002blcTldl ). 
IRAS 23151+5912 shows H2O maser emission and it is not 
detec t ed in the cm band at 8 GHz down to 1 mJy dTofani et al.1 
1995USridharan et alJl2002h . thus no significant hy percompact 



Hn region has formed yet ( Sridharan et al. 2002b . Based on 



single-dish 1.2 mm dust continuum observation s the total gas 
mass within th e region is estimated to ~6OOM ( Beuther et alj 
2002bl l2005ah . Observations of the CO(2-l) and SiO(2-l) 
lines with the IRAM 30 m telescope re veal a bipolar outflow 
with an entrained gas mass of ~ 20 M G dBeuther et alj|2002cb . 
The derived outflow rate of 1 x 10~ 3 M yr~' translates into an 
estimate of the accreti on rate of the same order of magnitude 
dBeuther et al.l l2002ch . Recent interferometric high-spatial- 
resolution SiO(2-l) observations reveal a quasi-parabolic blue- 
shifted outflow cone being con s istent with a wide-angle wind 
dOiu et al.ll2007b . IWeigelt et all d200fj|) speckle-imaged the re- 
gion in the near-infrared K-band, and the nebulosity they find 
coincides with the blue outflow wing. Their data rather suggest 
that the outflow features are produced by a precessing jet. In 
addition. lWeigelt et al. (2006) find an embedded K-band source 
IRS 1 which they propose to be the driving source of the molec- 
ular outflow. Furthermore, NH3 inversion lines from the (1,1) 
to the (4,4) transitions have been detected with the Effelsberg 
100 m telescope (unpublished data, Beuther priv. comm.), the 
excitation temperature of the (4,4) line is about 200 K. The NH3 
(1,1 )/(4,4) line ratio of only ~ 2 i ndicates warm and dense gas 
at the core center. Very recently, iGarav etaL ( 2007 ) detected 
an unresolved hypercompact Hn region at a sub-mJy flux level 
toward IRAS 23151+5912. Furthermore, they detected about 
30" apart a more evolved cometary Hn region. While the latter 
may contribute pa rt of the total luminoisty of the region (~20%, 
Garav et al.l 120071) . the main power-house of the region very 
likely remains the HMPO IRAS 23151+5912. Therefore, our 
target source comprises all features of a luminous High-Mass 
Protostellar Object at an early evolutionary stage. Assuming 
that the source is still in its a ccretion phase, it has not yet 
reached its final luminosity (e.g., Srid haran et alj2002b . and the 
total luminosity of the forming cluster will even exceed 10 5 L 
during its ongoing evolution. 

2. Observations 

We have observed the HMPO IRAS 23151+5912 with the 
Submillimeter Array during three nights between May and 



November 2005. We used three different array configurations 
(compact, extended, very extended) with unprojected baselines 
between 16 and 500 m, resulting at 875 /mi in a projected base- 
line range from 8 to 591 LI. The chosen phase center was R.A. 
[J2000.0]: 23 /, 17"'21. 5 and Deck [J2000.0] +59°28'48."995. 
The velocity of rest is v/ sr ~ -54.4 km s -1 . 

For bandpass calibration we used 3C279 and 3C454.3. The 
flux scale was derived in the very extended configuration from 
observations of 3c279. For the extended and compact configu- 
rations, we used 3C454.3 for the relative scaling between the 
various baselines and then scaled that absolutely via observa- 
tions of Uranus. The flux accuracy is estimated to be accu- 
rate within 20%. Phase and amplitude calibration was done 
via frequent observations of the quasars J0102+584 and BL 
Lac, about 13.5° and 20.7° from the phase center. The zenith 
opacity t(348GHz), estimated from the measurements of the 
NRAO tipping radiometer located at the Caltech Submillimeter 
Observatory, varied between the different observations nights 
between ~0.1 and ~0.3 (scaled from the 225 GHz measure- 
ment, r(348GHZ) ~ 2.5 x r(225GHz)). The receiver operated 
in a double-sideband mode with an IF band of 4-6 GHz so that 
the upper and lower sideband were separated by 10 GHz. The 
central frequencies of the upper and lower sideband were 348.2 
and 338.2 GHz, respectively. The correlator had a bandwidth 
of 2 GHz and the channel spacing was 0.8 125 MHz. Measured 
double-sideband system temperatures corrected to the top of 
the atmosphere were between 100 and 1000K, depending on 
the zenith opacity and the elevation of the source. Our sen- 
sitivity was dynamic-range limited by the side-lobes of the 
strongest emission peaks and thus varied between the line maps 
of different molecules and molecular transitions. This limita- 
tion was due to the incomplete sampling of short uv-spacings 
and the presence of extended structures. The lcr rms for the 
velocity-integrated molecular line maps (the velocity ranges for 
the integrations were chosen for each line separately depend- 
ing on the line-widths and intensities) ranged between 59 and 
80 mJy beam -1 (TableQ]). The average synthesized beam of the 
spectral line maps was 1" x 0.8" (TableQ]). The 875 /mi submm 
continuum image was created by averaging the line-free parts 
of both sidebands. The lcr rms of the submm continuum im- 
age was ~4.9 mJy/beam, and the achieved synthesized beam 
was 0.57" x 0.48" (PA. 86°). The initial flagging and calibra- 
tion was done with the IDL superset MIR originally developed 
for the Owens Valley Radio Observatory dScoville et al.| [l993) 
and adapted for the SMA0. T he imaging a nd data analysis were 
conducted in MIRIAD (Saul t et al.ll 1995b . 



1 The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics, and is funded by the 
Smithsonian Institution and the Academia Sinica. 



3. Results 

3.1. 875/um continuum emission in IRAS 231 51 +591 2 

Figure Q] shows the 875 /mi continuum image extracted from 
the line-free part of the data-cube. The spatial resolution of 
0.57" x 0.48" corresponds to -3000 AU at the given distance 
of ~5.7kpc. The peak and integrated fluxes of the main peak 
1 are about double of those from the minor peak 2 (Table 0. 



2 The MIR cookbook by Charlie Qi can be found at 
http://cfa-www.harvard.edu/~cqi/mircook.html 
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Table 1. Peak fluxes densities, rms and synthesized beams of 
integrated line images in Fig. [3] 



Line 


S peak 
mJy 
beam 


rms 

mJy 
beam 


Beam 
"(P. A.) 


CH 3 OH(7 3 ,5 - 6 2>4 ) 


1014 


79 


0.8x0.67(19) 


CH 3 OH(7 4 , 3 - 6 4 , 3 ),v, = 1 


617 


80 


1.0x0.81(19) 


C 34 S(7 - 6) 


411 


64 


1.08 x0.87(18) 


SiO(8 - 7) 


433 


59 


1.14x0.89(80) 


S0 2 (24 2 , 22 - 23 3 ,2i) 


826 


60 


1.07 x0.84(80) 


34 SO(8 8 - 7 7 ) 


977 


67 


1.08 x0.87(19) 



However, one of the H2O maser positions is clearly associated 
with submm peak 2 and none with submm peak 1. The sec- 
ond H2O maser position is right between both submm peak 
positions. Int e restin gly, the infrared source IRS1 detected by 
Weigelt et al ] d2006h is offset by about 0.9" north-east of the 
main submm peak position. The positional accuracy of the in- 
terferometric submm data is estimated to be ~ 0.1". Contrary 
to that, the astrometry of the infrared data is based on the 
2MASS astrometry, and this is estimated in this field to be not 
better than 1" (Preibisch, priv. comm.). Therefore, we cannot 
unambiguously decide whether IRS1 and the submm peak 1 
are the same or different sources. Based on the estimated H2 
column densities of the order 10 24 crrT 2 (corresponding to vi- 
sual extinctions of the order 1000, see below), one would not 
expect to detect near-infrared emission from the submm peaks. 
However, arguments based on the outflow emission may ex- 

20071) detected 



plain this (see 34.21) . Furthermore, IQiu et al. 
more extended 3.4 and 1.3 mm emission west of the peak. 
Although their 1 .3 mm feature is only at a 2 - 3<x level, the mor- 
phological similarity with the 3.4 mm data and the infrared/SiO 
outflow stimulates the speculation that this mm emission may 
be caused by the outflow. Extrapolating the 1 .3 mm flux density 
of ~ 4.25 mJy with a typical jet-frequency dependence of v 6 
to our continuum frequency of 875 /vm, the expected submm 
flux density would be S(875/mi) ~ 5.4 mJy, only barely above 
our lcr noise level. Therefore, it is no surprise that we cannot 
detect these features in the submm data. 

Table [2] lists the absolute source positions, their 875 //m 
peak fluxes and the integrated flux density approximately as- 
sociated with the two sub-sources. Assuming that the observed 
submm continuum emission is due to optically thin dust emis- 
sion, w e can calculate ga s masses and gas column densities fol - 
lowing iHildebrandl d 19831) and iBeuther et all d2002bi l2005ah . 
Since IRAS23151+5912 is in a very early evolutionary stage 
we use the dust tempe rature derived from the IRAS data (68 K, 
Sridharan et al.1 120021) and a dust opacity index ft = 2, cor- 
responding to a dust opacity per unit mass of /e(875yum) ~ 
0.8 cm 2 g _1 . Given the uncertainties in ft and T, we estimate the 
masses to be correct within a factor 2-5. Table [2] gives the de- 
rived masses and column densities separated for each of the two 
peaks. While the two core masses vary between 26 and 12 Mq, 
the column densities are of the order 10 24 cm 2 , corresponding 
to extremely high visua l extinctions o f the order 1000 (A,, = 
Nh/0.94 x 10 21 . iFrerking et alJfT982h . Integrating the whole 
submm continuum emission shown in Fig. [TJ(,Sj nt ~404mJy), 



Grey: 875,u.m continuum 

Triangles: H 2 maser (Beuthcr et al. 2002) 
IR sources (Weigelt ct al. 2006) 




20 5 



R.A. [J2000] 



Fig. 1. The grey-scale with contours shows the 875 fim con- 
tinuum emission. The contours go from +18 (3cr) to +98% 
(step +10%) from the peak flux of 80.8 mJy/beam. The two 
main submm peaks are marked and the triangles and stars 
show the posi t ions of H?Q maser and nea r -infrar ed peaks from 
Beuther et all d2002dh and IWeigelt et aD d2006h . The synthe- 
sized beam of 0.57" x 0.48" is shown at the bottom right. 

the total detected mass is ~38 M . Comparing the integrated 
flux density with the peak flux observed with the SCUBA ar- 
ray on the JCMT (~ 15" resolution) at 850/mi of 1890mJy 
dWilliams et al.ll2004 . approximately 80% of the total flux is 
filtered out in our interferometer observations^. The 3<x rms of 
the continuum image of ~14.7mJy corresponds to a mass sen- 
sitivity of ~ 1 .4 M at 68 K. Assuming lower temperatures off- 
set from the peak, the 3<x value corresponds to a higher mass 
sensitivity, e.g., -6.4 M at 20 K 

3.2. Submm spectral lines in IRAS231 51+5912 

The spectral bandwidths of 4 GHz in total revealed 27 spec- 
tral lines (Figf2]and Table [3). We detected a series of ground- 
state and vibrationally-torsionally excited CH3OH lines, some 
sulphur-bearing species ( 34 SO, SO2, C 34 S), the outflow- and 
shock-tracing SiO and a few more complex molecules like 
CH3CH2CN and CH3OCH3. The excitation temperatures of es- 
pecially the CH3OH lines cover a broad range from approxi- 
mately 40 to 360 K (TableEJ. 

Excluding the two spectral lines from the more complex 
molecules CH 3 CH 2 CN and CH3OCH3 that were too weak for 
imaging, we were able to image the molecular emission of 
all other species (Figf3]l. The three spectral line images from 
CH3OH, SO2 and 34 SO all peak in the near vicinity of the main 
submm peak and show emission approximately associated with 



3 Using only the compact configuration data, still approximately 
60% of the single-dish emission remain filtered out. 
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Table 2. Submm continuum source parameters 



Source 


R.A. 


Dec. 


S peak 


5 int 


M 


N 




[J2000] 


[J2000] 


r_mJy_i 

L beam * 


[mJy] 


[M G ] 


[10 24 cnT 3 ] 


Peak 1 


23:17:20.916 


59:28:47.65 


81 


282 


26 


2.3 


Peak 2 


23:17:20.800 


59:28:47.17 


44 


129 


12 


1.2 




Table 3. Line parameters 



348 

Frequency [GHz] 




Frequency [GHz] 

Fig. 2. Upper and lower sideband (top and bottom panel) spec- 
tra taken in the uv-domain on a short baseline of ~14m. The 
chosen spectral resolution is 2 km/s. 



Freq. 


Line 


E„ 


la- 


GHz 




K 


r mJ y i 

L beam 


337.348 


CH 3 CH 2 CN(383,36 - 37 3 , 35 ) 


328 




337.397 


C 34 S(7-6) 


65 


60 


337.582 


34 SO(8 8 - 7 7 ) 


86 


70 


337.642 


CH 3 OH(7 1 , 7 -6 1 , 6 )E(v,=l) 


356 


80 


337.722 


CH 3 OCH 3 (7 4 , 4 - 6 3 , 3 )EE 


48 




338.125 


CH 3 OH(7 ,7 - 6 ,s)E 


78 




338.306 


S0 2 (14 4 . 14 - 18345) 


197 




338.345 


CH 3 OH(7,, 7 - 6i, 6 )E 


71 




338.409 


CH 3 OH(7 , 7 - 6 , 6 )A 


65 




338.457 


CH 3 OH(7 5 , 2 - 6 5 ,i)E 


189 




338.486 


CH 3 OH(7 5 , 3 - 6 5 , 2 )A 
CH 3 OH(7 5 2 - 654 )A 


203 
203 




338.513 


CH 3 OH(7 4 4 - 64.3 )A" 
CH 3 OH(7 4)3 - 6 4 , 2 )A 
CH 3 OH(7 2 , 6 - 6 2 , 5 )A- 


145 
145 
103 




338.530 


CH 3 OH(7 4 , 3 - 6 4 , 2 )E 


161 




338.541 


CH 3 OH(7 3 . 5 - 6 3 4 )A + 


115 




338.543 


CH 3 OH(7 3 4 - 6 3 3 )A- 


115 




338.560 


CH 3 OH(7 3i5 - 6 3 , 4 )E 


128 




338.583 


CH 3 OH(7 3 , 4 - 6 3 , 3 )E 


113 




338.612 


SO 2 (20 U9 - 192,1,) 


199 




338.640 


CH 3 OH(7 2 . 5 - 6 2 , 4 )A 


103 




338.722 


CH 3 OH(7 2 . 5 - 6 2 , 4 )E 


87 


80 


338.723 


CH 3 OH(7 2 6 - 6 2 , 5 )E 


91 




347.331 


28 SiO(8-7) 


75 


60 


348.388 


S0 2 (24 2 , 22 - 23 3>21 ) 


293 


58 


349.107 


CH 3 OH(14 U3 - 14 , 14 ) 


43 





the region of submm continuum emission. This is partly dif- 
ferent for the three other line images. The integrated SiO(8- 
7) emission shows a relatively complex morphology which is 
largely caused by the molecular outflows in this region (see 
14.1b . The virbrationally-torsionally excited CH3OH emission 
shows a double-peaked integrated map around the main submm 
peak 1 and may be attributed to rotation from a larger-scale 
toroid around the most massive protostellar obect(s) (see !4.4l >. 
Furthermore, the C 34 S emission is centered offset from the 
submm continuum peaks and shows relatively strong more ex- 
tended emission reminiscent to the recent obser vations of the 
hot molecular core G29.96 dBeutheretalJl2007bh . 



Compared with hot cores like Orion -KL (IBeufher et al 



2005bl) or G29.96 dBeuther et alJ l2007bh IRAS 23151+5912 
exhibits relatively few spectral lines, however, it shows more 
lines than lower luminosity HMPOs like IRAS 05358+3543 
(Leurini et al. subm. to A&A). In a follow-up paper we present 
a more detailed comparison between all these sources with 
a special emphasis on the chemical diversity and a potential 
chemical evolutionary sequence (Beuther et al. subm. to A&A). 



4. Discussion 

4.1. Multiple molecular outflows 

The SiO(8-7) emission allows us to investigate the molecular 
outfiow(s) in this region. However, in comparison to the SiO(2- 
1) emission with a velocity spread from -70 to -48kms _1 
( Oiu et al. 20071) . the observable total velocity spread in the 
7 = 8-7 line is only from -60 to -51 km s~ 1 , re latively close to 
the ve locity of the ambient gas of -54.4 km s" 1 (iSridharan et al 
2002). The difference between the 7 = 2-1 and 7 = 8-7 



lines can likely be explained by the higher excitation temper- 
ature of the latter line. The energy levels above ground E u /k 
are 6 and 75 K for the lower and higher excited lines, respec- 
tively. With typical g as temperatures in molec ular outflows of 
the order 30 K (e.g.. ICabrit & Bertouj [l992b . the 7 - 8 - 7 
line is much harder to excite and hence more difficult to de- 
tect at all. Separating spectral core from line wing emission is 
not an easy task when no pronounced broad line wings are ob- 
servable, however, since SiO is believe to be produced almost 
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CH 3 0H(7 1>7 -6 1>6 )v t =l 









S0 2 (24 2>22 23 3i21 ) 



34 



S0(8 B -7 7 ) 





-2 
Aa ["] 







Fig.3. Integrated images in IRAS 23151+5912 of the molecules identified in Figure[2]that are strong enough for imaging. The 
species are labeld in each panel and the synthesized beams are shown at the bottom left of each panel. The grey-scale shows 
the submm continuum em ission as presented i n Fig. Q] Triangles and stars mark the H2O maser and near-infrared peaks from 
Beufher et al.1 d2002dl) and Weigelt et al. ( 2006h . The contour levels start at ±3<x and continue in +1.5<x steps, the lcr values are 
given in Table [3] 



exclu sively by sputtering from dust grains (e.g. JSchilke et al 
19971) . it is likely that most of the SiO emission is caused 



by the molcular outflow. Therefore, one can produce outflow 
maps from SiO data going closer to the ambient gas veloci- 
ties than it is feasible for more common molecules like CO. 
Hence, we imaged the blue- and red-shifted SiO(8-7) emis- 
sion in the intervals [-60,-56] km/s (blue) and [-54,-51] km/s 
(red), respectively. Figure |4] presents this SiO line wing maps 
as overlay on the submm continuum emission. The blue- and 
red-shifted emission are spatially distinct, but the morphology 
is more ambiguous than would be expected fro m the previous 



mm and near-infrared outflow observation s (e.g jBeufher et al 



l2002d:IWeigelt et al.l2006llOiu et al.l2007l) . These previous ob- 
servations clearly identify one outflow in approximately east- 
west direction emanating from the region containing the main 
submm peak and the infrared source IRS 1 . This east- west out- 
flow can be recognized in the new SiO(8-7) data as well, al- 
though one should note that the geometric center of the blue- 
and red-shifted peak is not exactly the main submm contin- 
uum peak but offset about 0.5" to the west. Furthermore, there 
is considerable additional SiO blue and red-shifted emission 



north and south of the continuum sources, respectively. This 
emission is offset from the dust continuum peaks and hence 
unlikley be caused by the ambient gas. Since we have at least 
two embedded protostellar sources, we suggest that we are 
witnessing two outflows emanating from th e region as well . 
Reinspecting the SiO(2-l) observations from lOiu et al. ( 2007 ). 
one finds that their lower excitation data show an additional 
red-shifted isolated peak toward the south as well. Our SMA 
data were observed in three different configurations having a 
relatively good uv-sampling and hence s patial coverage, how- 
ever, the southern peak from lOiu et alj ( 12007b is close to the 



edge of our primary beam and thus difficu lt to detect any way. 
It is likely that the SiO observations by lOiu et al.1 d2007l) fil- 
tered out more of the extended emission which could cause the 
relative isolation of the southern SiO peak in their data. 



4.2. The outflow driving sources 

Based on a fa i nt jet-l ike infrared feature emanating from IRS 1 , 
Weigelt et al. (2006) propose that this should be the driv- 
ing source of the east-west molecular outflow. In contrast to 
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Fig. 4. The full and dashed contours show the blue- and red- 
shifted SiO(8-7) emission for IRAS 23151+5912 integrated 
from -60 to -56 and from -54 to -5 1 km/s, respectively. Thick 
and thin lines present positive and negative features, the latter 
ones being due to the incomplete uv-coverage of interferome- 
ter datasets. The two lines highlight the two tentatively identi- 
fied molecular outflows, and the synthesized beam is shown at 
the bottom-right. The grey-scale shows the submm continuum 
emission as presented in Fig. [T] Triang les and stars mark the 
H? Q maser and near-infr ared peaks from lBeuther et al.l |2002d) 
and lWeigelt et al.l(12006l) . The contour levels go from 18 to 98% 
(step 10%) of the peak emission with blue- and red-shifted peak 
values of 552 and 445 mJy beam -1 , respectively. 



duced a position-velocity diagram through submm peak 2 at a 
position angle of 160 degrees from north, proximately the ori- 
entation of the 2nd proposed outflow (Pig. [5j. While it is less 
clear for the blue-shifted emission, the red-shifted side shows 
a clear increase of velocity with distance from submm peak 
2, resembling the Hubble-law for molecular outflows (e.g., 
Rich er et al.l l2000). Based on the morphology and the position- 
velocity structure of the north-east south-west outflow, we sug- 
gest that this outflow is driven by the embedded source within 
submm peak 2. 



° ( 



-1 .90 ) PA= 




-56 -54 
velocity (km/s) 

Fig. 5. SiO(8-7) position-velocity digram of the north-east 
south-west outflow centered on submm peak2 (offset -1.5"/- 
1.9" from the phase center) with a position angle of 160 de- 
grees from north. Negative offsets are going north-west. The 
contour levels go from +0.122 to +0.732 Jy beam -1 in steps of 
+0.122Jybeairr'. 



that, one would expect that powerful outflows like that in 
IRAS 23151+5912 are driven by younger and more deeply em- 
bedded sources like the submm peak 1 . While we cannot astro- 
metrically decide whether the two sources are the same or not, 
the outflow associations indicate that the infrared and submm 
emission may indeed emanate from the same massive protostel- 
lar object. Furthermore, due to scattering effects within embed- 
ded objects it is possible that the emission peak shifts at near- 
infrared wavelengths in the direction of the stee pest density 



5ep_e 

gradient which is towa rd the outflow cavity (e.g., Eisner et all maser emission from low-mass regions dWilking et alj 1 1994 



The second H2O maser position right between the two 
submm peaks may either be caused solely by shock interac- 
tions of a molecular outflow with the ambient gas, or it could 
be associated with a low-mass source which is below our de- 
tection threshhold (3cr continuum sensitivity corresponds to 
~1.4M at 68 K or a higher value of ~6.4M at lower tem- 
peratures of 20 K, jj3.U and hence not possible to be identified 
as an additional weak submm peak. Although H2O masers are 
more often found towar d massive star-forming regions (e.g., 
Forste r & Caswell 119991) . there exist seve ral examples of H7O 



20051 IStarketail20 06). Therefore, while an outflow cavity al 
lows to penetrate deeper into the core than possible in spherical 
symmetry, the scattering effects could cause even a real offset 
between the infrared and submm source. Hence, although the 
near-infared source IRS 1 and the submm peak 1 may be two 
separate sources, the data are also consistent with both emis- 
sion peaks potentially being due to the same embedded young 
and massive protostellar source that could drive the east-west 
outflow. 

The association of one of the maser peaks with submm 
peak 2 potentially indicates that this source could be driving 
a molecular outflow as well. To check this hypothesis, we pro- 



Claussen et al.lll998l) as well as from lower-mass companions 



in regions of high-mass star formation dHunter et al . 1999). 



4.3. The CH 3 OH line forest 

The detection of 19 CH3OH lines with excitation tempera- 
tures between 43 and 356 K allows the temperature to be es- 
timated toward the two submm continuum peaks. We modeled 
the whole lower sideband CH3OH spectrum toward both po- 
sition in the local thermodynamic equilibrium approximation 
using the XCLASS superset to the CLASS software developed 
by Peter Schilke (priv. comm.). This software package uses the 
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line catalogs from JPL and CDMS dPovnter & Picked 1 1985 
Mull er et alJliooTh . The important source dependent model in- 
put parameters are the temperature, the column density, the 
source size and the line-width. While the latter one is relatively 
easy to determine from the data, the former three parameters 
are partly degenerate and various permutations give reasonable 
fits. Figure [6] shows the observed lower sideband spectra with 
reasonable model spectra overlayed, the exact model input pa- 
rameters are stated in each panel. The most important value 
for the characterization of the two cores is the temperature. 
Therefore, we went through the parameter space of mainly tem- 
peratures and column density to estimate the uncertainty. The 
derived temperatures with approximate error-bars for the two 
positions are T(Peakl) ~ 150+50 K and 7\Peak2) ~ 80+30 K. 
Although the error-bars are large, as expected from the detec- 
tion of the vibrationally-torsionally excited line at the location 
of submm peak 1, the data indicate that temperatures toward 
that peak are likely to be higher. 



flows. Similar to so me other recently in vestigate d regions (e.g.. 



IRAS 18089-1 732 iBeufher et al.lEo05d G29.96 lBeuther et al 




3.38 10 3.385 
Rest Frequency (MHz) 

Fig. 6. Lower sideband spectra toward the two submm peaks. 
The grey-scale and dotted contours show the orginal spec- 
trum at 1.1" x 0.9" resolution, and the thick solid lines present 
synthetic spectra assuming Local Thermodynamic equilibrium 
with the parameters stated in the figure. 



4.4. Rotational signatures 

One of the missing links in understanding high-mass star for- 
mation is still t he detailed understan d ing of massive accre - 



tion disks (e.g., ICesaroni et al.l 12007c iBeuther et al. 2007a) 



2007b), most of the lines do not show any coherent veloc- 
ity signature indicative of rotation. The only exception is the 
vibrationally-torsionally excited methanol line CH30H(7ij - 
6i,6)vr = 1 which shows a blue-red velocity shift in north-south 
direction around submm peak 1 (Fig.|7j. Since CH3OH v, = 1 
has an upper level excitation temperature of ~356 K (TableO it 
traces only the warm gas likely excited by the central growing 
massive protostar. Furthermore, this velocity gradient is per- 
pendicular to the main outflow ( {34.11 & Fig. [U, which likely 
indicates rotational signatures caused by an underlyin g collaps- 
ing rotational core-disk system (e.g JCesar oni 2005). The sep- 
aration between the two blue- and red-shifted CH3OH v, = 1 
peaks is approximately 4300 AU. In addition to the relatively 
large size of the rotating structure, the moment 1 map and 
position-velocity digram in Figure [7] do not show a Keplerian 
velocity pattern as observed for many low-mass sources (e.g., 
Simon et alj|2000l) . A Keplerian profile would show the largest 
velocities close to the center with a velocity decrease propor- 
tional to r whereas here we see even a velocity increase 
with distance from the center. Therefore, this structure is not a 
Keplerian accretion disk but it rather resembles a larger-scale 
rotating toroid that may harb or a genuine ac cretion disk at its 
so far unresolved center (e.g.. lCesaronill2005l) . 

Although this structure is apparently not in Keplerian ro- 
tation, one can nevertheless estimate an approximate binding 
mass of this rotating structure M mt assuming equilibrium be- 
tween the rotational and gravitational forces at the outer radius. 
Then one gets 



2 
v r 



M rn , = 



M rot [M ] = 1.13 10~ 3 x v 2 [km/s] x r[AU] 



(1) 

(2) 



where r is the radius of the rotating structure (~2150 AU, half 
the peak separation in Fig|7]top panel), and v the Half Width 
Zero Intensity (HWZI) of the spectral line (~3 kms~', half the 
velocity regime used for the integration in Figj7] top panel). 
EquationsQ]&[2]have to be divided by sin 2 (/) where i is the un- 
known inclination angle between the plane of the rotating struc- 
ture and the plane of the sky (i = 90° for an edge-on system). 
With the given values we get a mass M ml of ~ 22/(sin 2 (/))M , 
similar to the mass value derived for Peak 1 from the submm 
continuum emission (Table [2]). Assuming an inclination an- 
gle of 45 degrees, the rotationally supported mass would be 
even 44 M . Although higher than the gas mass derived from 
the dust continuum emission, this is still within the error mar- 
gins of the dust continuum calculation which is usually esti - 



mated to be correct within factors 2-5 ( Beufher et alJl2002bl) . 



Therefore, we searched the spectral-line data-cube for rota- 
tional signatures perpendicular to any of the molecular out- 



Assuming that the submm continuum emission largely stems 
from the disk and/or rotating envelope, this implies that the 
disk/envelope mass has to be of the same order than the mass of 
the central embedded protostar. This is different to typical low- 
mass protostar-disk systems where the disk-mass is usually 
negligable to the protostellar mass. However, one has to note 
that such typical low-mass systems are usually more evolved 
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and have already dispersed much of their envelopes whereas 
IRAS 23151+5912 is still deeply embedded within its natal 
core. In addition to this, as shown above, the rotating structure 
is not in real Kelplerian rotation, and therefore the assumption 
of equilibrium between gravitational and rotational force does 
exactly hold from the beginning. 

5. Conclusion and Summary 

We present submm line and continuum observations ob- 
tained at sub-arcsecond spatial resolution with the SMA to- 
ward the very young 10 5 L High-Mass Protostellar Object 
IRAS 23151+5912. The submm continuum emission is re- 
solved into two sub-sources, each of them likely driving one of 
the molecular outflows observed in SiO(8-7). The main submm 
continuum peak is approximately 1" offset from a recently 
identified near-infrared sour ce which is suggeste d to drive one 



of the molecular outflows (Wei gelt et al.l 120061) . With these 



data, we cannot differentiate whether this offset is real indicat- 
ing the presence of an additional source, or whether it is an ob- 
servational artifact caused by the relatively poorer near-infrared 
astrometry. The detection of an intermediate dense spectral line 
forest with 27 lines within the 4 GHz bandpass allows the in- 
vestigation of additional physical and chemical properties. The 
19 CH3OH lines provide temperature estimates toward the two 
main submm continuum peaks to 7\Peakl) ~ 150 ± 50 K and 
7\Peak2) ~ 80 ± 30 K, respectively. Searching the dataset for 
velocity signatures indicative of rotation around the central pro- 
tostar, only a vibrationally-torsionally excited CH3OH line ex- 
hibits a velocity gradient perpendicular to the main molecular 
outflow. With a size-scale of a few 1000 AU and no Keplerian 
velocity profile, we are likely not observing a real accretion 
disk but rather a larger-scale rotating toroid which may harbor 
a genuine accretion disk at its very center. 

Setting IRAS 23151+5912 into context with lower lumi- 
nosity HMPOs and higher luminosity UCHns, its observed 
physical and chemical parameters (e.g., line forest, outflow 
morphology, rotation signatures) are a mixture of what is ob- 
served in the two other classes. This supports the idea that this 
region is in a kind of transition phase between the HMPO and 
the UCHn region stage. In a follow-up paper, we will present 
a comparison of the chemical properties of this regions with a 
few sources representative of different evolutionary stages and 
luminosities observed all in the same spectral setup with the 
SMA (Beuther et al. subm. to A&A). 
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i ( -0.60 , -1.30 ) PA= 0.00 




Fig. 7. CH 3 OH(7i ;7 - 6i,6)v, = 1 emission for 
IRAS 23151+5912. Top-panel: The full/dashed contours 
show the blue-/red-shifted emission integrated from -58.5 to 
-55.5 and from -54 to -52km/s, respectively. Thick and thin 
lines present positive and negative features, the latter ones 
being due to the incomplete uv-coverage of interferometer 
datasets. The grey-scale shows the submm continuum emission 
as presented in Fig.Q] Tri angles and stars mark t he H?Q maser 



and n ear-infrared peaks (Be uther et al.l 2002d; Weig elt et al 



120061) . The contour levels start at the +2cr levels and continue 
in +lcr steps with the lcr levels of 87 and lOOmJybeairT 1 for 
the blue- and red-shifted maps, respectively. Middle panel: 
The grey-scale shows the 1st moment map and the contours 
outline the submm continuum emission as presented in Fig.Q] 
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